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Abstract 

In this paper we obtain a heavy meson bound state equation from the heavy 
quark equation of motion in heavy quark effective theory (HQET) and the 
heavy meson effective field theory we developed very recently. The bound 
state equation is a covariant extention of the light-front bound state equation 
for heavy mesons derived from light-front QCD and HQET. We determine the 
covariant heavy meson wave function variationally by minimizing the binding 
energy A. Subsequently the other basic HQET parameters Ai and A2, and 
the heavy quark masses nii, and rric can also be consistently determined. 
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Typeset using REVT^ 



I. INTRODUCTION 



Very recently, we have developed an effective field theory to describe the nonperturba- 
tive QCD dynamics of heavy mesons [0,0 . This theory has incorporated in it heavy quark 
effective theory (HQET), and is thus consistent with heavy quark symmetry (HQS) in the 
infinite quark mass (mg —>■ oo) limit. In this effective field theory, a heavy meson is con- 
sidered as a composite particle consisting of a reduced heavy quark (heavy quark in the 
ruQ — >■ oo limit) coupled with the light degrees of freedom. A structure function \I^(f -pq), 
corresponding to the wave function of the heavy meson bound state, is introduced to de- 
scribe the quark-meson vertex in the effective Lagrangian. By combining such an effective 
Lagrangian for the composite heavy mesons with the l/mq expansion of the heavy quark 
QCD Lagrangian, we can systematically evaluate various l/mq corrections to heavy meson 
properties in a covariant field-theoretic framework [Q. 

In our previous works |l|J^, the heavy meson wave function \l/(f ■ Pg) is taken to be a 
phenomenological function. Specifically it is assumed to be a Lorentzian function, \E'(f -p^) = 
1 / (v-pq+uj)"', in which the parameters u and n are determined phenomenologically by fitting 
to the B meson decay constant fs', hence its relation to nonperturbative low-energy QCD 
dynamics is not clear. In this paper, we attempt to determine the heavy meson wave function 
from a QCD-based bound state equation. The bound state equation of heavy mesons is 
obtain from the heavy quark equation of motion in HQET which is incorporated in our 
effective field theory. This heavy meson bound state equation is a covariant extension of the 
light-front bound state equation we obtained recently from light-front QCD and HQET [^. 
With the bound state equation for \I'(f -pg), we can use a variational procedure to determine 
its soluton. 

The paper is organized as follows: In Sec. II, we derive from the heavy quark equation of 
motion in the heavy quark limit and our effective heavy meson field theory the heavy meson 
bound state equation. Then a variational method is used to determine the heavy meson wave 
function \l/(f ■ pg) by minimizing the reduced mass (binding energy) of the heavy meson A. 
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In the literature, A has only been estimated either from heavy meson mass spectrum or from 
inclusive heavy meson decay data A direct calculation of A has not been carried out to 
date. In Sec. Ill, we consistently calculate the other two basic HQET parameters Ai and 
A2 directly from the operators Oi and O2 in HQET. Having self-consistently calculated Ai 
and A2, we can then extract from experimental heavy meson masses the heavy quark masses 
[nib, rric) at the low energy scale where our effective field theory applies. Such a procedure 
establishes a link between our heavy meson effective field theory and nonperturbative QCD 
dynamics. Finally a summary is given in Sec. IV. 

II. HEAVY MESON BOUND STATE EQUATION IN EFFECTIVE FIELD 

THEORY 

A. Heavy meson effective field theory 

Before deriving the heavy meson bound state equation, we shall briefly review the heavy 
meson effecive field theory we developed recently In heavy quark effective theory 

(HQET) the QCD Lagrangian is expanded in powers of l/mQ-. 

00 / 1 \n_ 

+ E — K{^:^±){-^v ■ Dr-\i]^^)K 

„=1 K^^Q/ 

= C^ + (2.1) 

where the reduced heavy quark field is defined from the full heavy quark field Q by 

g(a;) = e-^™'3^-^/i, + 0(l/mQ) , ^^h, = K , (2.2) 

is the heavy quark four velocity which is defined to be equal to that of the parent heavy 
meson, and = D'^ — v^v ■ D. It is easily seen that the leading term Cq = h^iv ■ Dh^ 
in Eq.(|2.1|) possesses a SUf{Nf) ® SUs{2) flavor-spin symmetry (or simply heavy quark 
symmetry, HQS), whereas the symmetry breaking part Cmg contains l/mq corrections to 
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The effective Lagrangian we have constructed to describe the low energy dynamics of 
heavy mesons reads 



Ces = Kiv ■ dhjj + q{ip - mq)q - -F^Ta^^u 



1 

4 

-(hJ'y^q^P^ - K-f^q^V^^ + h.c. 



OO / -1 \ " 



n=l V^qJ 

C^f + , (2.3) 



where Py and represent respectively the composite pseudoscalar and vector heavy meson 
fields which appear only as external states, 

A= lim (MH-mQ) (2.4) 

is their common residual mass in the heavy quark limit, 

qy = G^ivpg)q (2.5) 

represents collectively the light degrees of freedom in a heavy meson, where \I'(f -p) is heavy 
meson wave function mentioned earlier, pq is the light antiquark momentum, and G is the 
normalization constant given by 

J (27r)4 (A-vpg + lef (p^ - ml + te) ' ^ ' ^ 

In Eq. (|27^) , obeys HQS and defines the nonperturbative effective couphng of heavy 

mesons with their constituent heavy and light quarks, while Cmg, which contains all the 
l/rriQ corrections, is considered as a perturbation to C^^^ ■ A practical calculational scheme 
has been developed in terms of the Feynman diagrammatic rules summarized below [|T],3: 

(i) A heavy meson bound state {j>h = At>) couples to its constituent heavy and light 
quarks via a three-point vertex 



-iG^>{vpq)TH. (2.7) 



iG^!*{vp,)Th. (2.8) 

where Vp = 275 for pseudoscalar (if = P), Fy = — / for vector mesons {H = V), and Pq 
is the momentum of the hght degrees of freedom. Each vertex also carries a momentum 
conservation factor of (27r)^5^(Af — k — pg). 

(ii) The internal heavy quark and light quark propagators are, 

■ (2.9) 



2(vk + ie 



r :*'T.\ . (2.10) 



-Pq p\-m^^Ar le 

respectively, where k is the residual momentum of the heavy quark, and the constituent 
mass of the light quark. 

(iii) An integration factor / d!^p/ {2ttY is associated with each internal momentum p, and 
a factor (—1) should be included for each fermion loop. 

(iv) For all other lines and vertices not attached to heavy mesons (mostly from Cmg), 
standard Feynman rules apply. Hence l/mq corrections to various quantities can be calcu- 
lated using standard field theoretic perturbative approach 0. 

B. Heavy meson bound state equation 

In HQET, the reduced heavy quark field obeys the following equation of motion 

iv-Dh^ = 0, (2.11) 

where = — igA^ is the QCD covariant derivative. 

To ensure that our effective field theory, Eq.( p.3|) , is consistent with low energy QCD 
dynamics of HQET, Eq.( p.l|) , we require that in the heavy quark limit, the heavy meson 
state \H{y)) satisfies 

{H^lKiv DK\H^) = 0. (2.12) 



Strickly speaking, this matrix element must be calculated in a nonperturbative manner of 
QCD, which is however beyond our present day knowledge of QCD. To circumvent this 
predicament, we shall adopt an effective gluon propagator which contains a long-distance 
contribution as well as the familiar perturbative one 

g^W'iq) = -ig^^5a, + i^7,'Al5\q)\ (2.13) 

where Qs is the QCD coupling constant, and the coefficient characterizes the confining 
strength. Note that the regularized infrared singular term 5'^{q) can be effectively considered 
as a regularization of the linearly confining interaction 1/q^ 0. Representing this effective 
gluon propagator by a double wavy line, we can depict the left hand side of Eq. ( ^.12[ ) 
pictorially by Feynman diagrams shown in Fig. 1. 

As discussed in [|T],0], from the constraint that it forbids on-mass-shell decay of the heavy 
meson into Qq, a possible form of the heavy meson wave function \l'(f ■ pq) is given by 

^>{vPq) = (K-VPq)^{vPq), (2.14) 

where the function ipiy ■ Pg) is regular at v ■ Pg = A. Then, applying the diagrammatic rules 
stated in Sec. IIA to Fig. 1, we can easily translate Eq.( |2.12| ) into the following expression. 



X 



{V ■Pg+ mg){V ■p'g + mg)+V ■ PgV ' p'g - P g ' p'g 



(p2 _m2)(p'J -m2) 
= 0, (2.15) 

where the color factor Cf = [N^ — l)/2Nc = 4/3 for Nc = 3. The above equation can be 
rewritten as an equivalent bound state equation for heavy mesons: 



(A - - pMv ■P.) = CfjT^A - ^-'AyiP. - P[ 



(2TlY\{Pg-p'^Y 
{V-Pg + rUg) {V ■ P'q +mq)+V- PgV -Pg-Pq-Pg 

{v ■ Pg + mg){p'l - ml) 
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ifiv-p'), (2.16) 



which is indeed a covariant generahzation of the hght-front bound state equation for heavy 
mesons we have previously derived 0], and it is also the heavy quark limit of the Bethe- 
Salpeter equation for heavy mesons. 

Using the normalization condition of ■ Pq), we can simplified Eq.( |2.15| ) to 

{pl - ml){p'l - ml) 

The physical meaning of the above result is clear. The heavy meson binding energy A 
defined by Eq.( p.4|) consists of two parts: One is the light quark kinetic energy contribution 
(f ■ Pq), and the other is the interaction energy of the heavy-light system. Note that, due 
to the effective gluon propagator (|2.13|) used, the interaction energy contains a confining 
contribution, as well as the color Coulomb contribution. 

C. Variational determination of heavy meson wave function {p{v ■ pq) 

In this section we shall determine the heavy meson wave function ip{v ■ Pq) through a 
variational procedure. The first step is to choose an appropriate trial wave function. As 
discussed in we demand that the wave function should be relatively well behaved, by 
which we mean that, when continued into the complex Pg plane, (i) (f{v ■ pg) is analytic 
except for isolated singularities, and (ii) it vanishes "fast enough" as \v- q\ oo. These two 
constraints also allow us to easily derive various relativisitic and/or light-front quark models 
results as special cases in our effective field theory 

In various existing light-front or relativistic quark models, two types of distribution 
functions have been widely adopted for phenomenological hadronic wave functions. One 
is the exponential(Gaussian)-type and the other the Lorentzian-type. The former can be 
excluded since it does not simultaneously satisfy the two constraints stated above On 
the other hand, we find that Lorentzian-type functions of the form, 
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Vn{v ■ p,) = — , (2.18) 

[V ■ Pg + uj — «e)" 

do meet our requirements [0,0, and will be taken as the working trial heavy meson wave 
function in this work. It is interesting to note that such a choice of (f{v ■ Pg) corresponds to 
a covariant light-front wave function 0] of the form 



»'-(-P,)-A..f^ ">4, (2.19) 

where Afi is a normalization constant, Pg = m^, and n > 4 is required to ensure convergence. 

For the sake of demonstration, we first consider the naive case of ignoring of the confining 
interaction between the heavy quark and the light degrees of freedom, i.e. = 0, so that 
color Coulomb interaction is solely responsible for binding the quarks together. We find 
that, for n = Q, a minimum with A = 0.192 GeV can be found at uj = 0.03. However a wave 
function with such a small u is unrealistic because it contains too little high momentum 
components. In other words, the wave function is too extended in configuration space, as 
would be expected since we have ignored the confining part of the inter-quark interaction. 

With an appropriate non-zero Ac, physically sensible results can be obtained. In what 
follows, we will take as = 0.44 and rriq = 0.22 GeV, which are close to the values used 
in our previous work p|. With the confinement strength taken to be A^ = 0.665 GeV, we 
determine the wave function parameter u by minimizing A. The results are listed in the 
Table I, where 

A = Eke + £'conf + ^coui = 0.426 ~ 0.430 GeV (2.20) 

which is not sensitive to the choice of n. This value of A is somewhat larger than our 
previous result of A = 0.33 ~ 0.34 GeV obtained in [^. The difference is mainly due to the 
confinement effect which is not considered there. 
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Table I. The wave function parameter lo nonperturbatively determined from bound state equation. 





LO 


{v-pi) 


-^conf 


-^coul 


A 


n = 6 


0.80 


0.443 


0.147 


-0.160 


0.430 


n = 7 


1.01 


0.440 


0.147 


-0.158 


0.427 


n = 8 


1.24 


0.438 


0.144 


-0.157 


0.426 



It is worth emphasizing that, in our previous work the parameter uj is fixed phe- 
nomenologically by fitting to fs = 180 MeV from lattice gauge calculations (it is experi- 
mentally unknown). In the present work, we have dynamically determine the wave function 
ipiy ■ p) by minimizing the blinding energy A from the QCD-based bound state equation. 
This is a significant improvement which helps to raise the predictive power of our heavy 
meson effective theory. 

III. PSEUDSCALAR AND VECTOR HEAVY MESON MASSES 

In this section, we shall use the bound state solution obtained in the previous section 
to calculate heavy meson masses up to order I/itlq. This will in turn yield more consistent 
results for Ai and A2 than those obtained in [Q. 

In the heavy quark limit, the heavy meson masses can be written as Mh = ttiq + A. The 
correction to Mh in HQET comes mainly from the leading HQS-breaking 1 / mg corrections 
[see Eq. Q]: 

£1 = h^D^fK + ^-^Ka,M^'K = 0^ + 02, (3.1) 

where a^i, = ^['y^,'ju] and G'^'^ = y[D^^,D^]. With these l/mg corrections included, the 
heavy meson masses can be expressed as 

= mQ + A- -^(Ai + dj^Aa), (3.2) 

where Ai comes from Oi, (1^X2 from O2 (see Fig. 3), and du = 3(— 1) for pseudoscalar 
(vector) mesons 0. In Eq. (|3.2| ), Ai parametrizes the common mass shift, and A2 accounts 
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for the so called hyperfine mass splitting between the vector and pseudoscalar heavy meson. 
They can be directly calculated once the heavy meson wave function ip[v ■ p) is determined 
from the variational approach presented in Sec. IIC. 

Expressions for Ai and A2 are essentially the same as those given in 0, except that here 
we use the effective gluon propagator ( ^.131 ), which is more appropriate for physics in the 
low energy regime. Thus we have 



4^/ 

9'' {^2 _ ^2\(w2 



X [{Vpq + mq)[{pq + p'g) ' p'g - V ■ {Pq + p'g)v ■ p'g] 

+ {V-P'g+ TUq) [{Pq +P'q)-Pq-V- {Pq + p'g)v ■ Pq]} , (3.3) 



and 



A, 



-CfG / - — TIT — hVyv-Va)^ [v-Vn 



(27r)4(27r)4 <? (p2 _ ^2);^^/^ _ ^, 

^^"^ +^A7^'Kl5\pq-p' 



{{Pq-Pq)' 

X [{Vpq+ mq)[{pq - p'q) ' p'q - V ■ {Pq - p'g)v ■ p'^] 

-{v-p'g + rrig) [{p'g -p'q)-Pq-v- {pg - p[)v -Pg]}. (3.4) 

Note that since the chromo- magnetic interaction vanishes at q = {pq = p'^), A2 is actually 
not affected by the confinement part of the gluon propagator. Now, we can use a self 
consistent method to determine uniquely the parameters of our effective theory (Ac, cj), and 
also those of HQET (A, Ai, A2, and mg). The procedure goes as follows. For a given Ac, 
we first determine the wave function parameter u by minimizing the binding energy A via 
Eq.( p.l7|) . Then we can evaluate Ai and A2 using Eqs. (|3.3| ) and ( p.4|) . From Eq. (p.2|) , the 
hyperfine mass splitting is given by 

2A 

AMvp = Mv-Mp = —, (3.5) 
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which, together with the experimental data AMb*b = 45.7 ± 0.4 MeV and AMo*d = 
142.12 ± 0.07 MeV |T0[, uniquely determines the heavy quark masses rrib and rric. Finally, 
knowing Mb = 5.279 GeV, we can recalculate A from Eq. (|3.2| ). This resulting A should be 
equal to that determined earlier from variation. If not, the process is repeated with a new Ac 
until we reach consistency. Through this self consistent procedure, we can for the first time 
determine uniquely all the basic HQET parameters m^, m^. A, Ai and A2 from QCD-based 
dynamics in the low energy regime. The final results are listed in the Table II. 



Table II. The wave function parameter uj is variationally determined by minimizing A, and all the other 
HQET parameters are calculated self consistently in the heavy meson effective theory with niq — 0.22 GeV 
and Us = 0.44. Ac, uj, nib, 'rric, and A are in GeV, and Ai and A2 in GeV^. 



"fniv -Pq) 


Ac 


UJ 


a: 


-Ai 


A2 




TOfc 


rric 


n = 6 


0.665 


0.805 


0.430 


0.228 


0.1116 


0.0459 


4.860 


1.570 


n = 7 


0.685 


1.070 


0.436 


0.247 


0.1114 


0.0459 


4.852 


1.568 


n = 8 


0.705 


1.336 


0.443 


0.264 


0.1112 


0.0459 


4.845 


1.556 


exp. 






0.39 ±0.11 ^ 


0.19 ±0.1 


0.1117* 


0.0457 ±0.004 


4.89 ± 0.05** 


1.59 ±0.02** 



*In the literature, the so-called "experimental" value A2 — 0.12 GeV^ is obtained by assuming nib = Mb- 



One would get A2 ~ 0.1117 GeV^ if the pole-mass to^"''' = 4.89 GeV is used in Eq. Q for AMb*b = 45.7 



MeV. 

** There are no experimental data for the constituent quark masses nib and TOc, and the pole masses are 
displaced |l3]. 



IV. SUMMARY AND CONCLUSION 

We know that HQS and confinement mechanism are crucial ingredients of nonpertur- 
bative QCD dynamics in heavy hadrons. However they are rarely incorporated in quark 
models in a fully covariant manner. In this work, we attempt to establish a link between 
quark model and low-energy QCD dynamics in an effective field theoretic approach, in which 
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the heavy meson wave function is determined variationally. Our effective theory preserves 
the simphcity of a conventional quark model; in fact, for heavy quark transitions, we can 
reproduce the covariant light-front quark model ||^ as a special case. The light-front quark 
model has been widely used to describe heavy meson structures. However, because its rela- 
tion to QCD is not clear, it is incapable of describing the hadronic mass spectrum. Moreover, 
the model is not fully covariant, so that it can only describe physical processes with space- 
like momentum transfers . In the present formalism, we can study hadronic bound state 
physics in a covariant framework; moreover, a connection to low-energy QCD dynamics has 
also been established. Thus our heavy meson effective theory, together with the nonpertur- 
batively determinated heavy meson wave function, provides a useful QCD based framework 
in which to study heavy quark physics in the low energy regime. 

We have adopted a self-consistent procedure to determine the heavy meson wave function, 
the binding energy A, as well as its l/mg corrections Ai and A2. /,From these HQET 
parameters, we can determine the heavy quark masses and in the low energy regime 
by fitting to the B* — B and D* — D hyperfine mass differences. For these static properties, 
our results are consistent with other studies on inclusive B-meson decays [^]. 

A great advantage of our formalism is that we can unambiguously study hadronic tran- 
sition processes at arbitrary momentum transfers. Using the heavy meson wave function 
determined in this paper, we can further calculate dynamically the heavy meson decay con- 
stant (which is usually considered as an input in various quark models for heavy mesons) 
and other hadronic properties such as the Isgur-wise function, DD*tt and DD*'~f coupling 
constants, and also their l/mq corrections. The results will be presented in forthcoming 
papers. 
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Figure Captions 

Fig. 1 Diagram for heavy meson bound states: (a) kinetic energy contribution and (b) 
effective one-gluon exchange contribution which contains the Columbus interaction 
plus confining interaction in heavy quark symmetry. 
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